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USING CONTINUUM MODELS, PART I 

A. V. Balakrishnan 


Abstract 


This paper is Part I of a two part report and represents a continuing systematic 
attempt to explore the use of continuum models — in contrast to the Finite Element 
Models currently universally in use — to develop feedback control laws for stability 
enhancement of structures, particularly large structures, for deployment in space. We shall 
show that for the control objective, continuum models do offer unique advantages. 

It must be admitted of course that developing continuum models for arbitrary struc- 
tures is no easy task. In this paper we take advantage of the special nature of current 
Large Space Structures - typified by the NASA-LaRC Evolutionary Model which will 
be our main concern — which consists of interconnected orthogonal lattice trusses each 
with identical bays. Using an equivalent one-dimensional Timoshenko beam model, we 
develop an almost complete continuum model for the Evolutionary structure. We do this 
in stages, beginning only with the main bus as flexible and then going on to make all the 
appendages also flexible — except only for the antenna structure. 

Based on these models we proceed to develop formulas for mode frequencies and 
shapes. These are shown to be the roots of the determinant of a matrix of small 
dimension compared with mode calculations using Finite Element Models, even though 
the matrix involves transcendental functions. The formulas allow us to study asymptotic 
properties of the modes and how they evolve as we increase the number of bodies whtch 
are treated as flexible - as we shall see the asymptotics in fact become simpler. 
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MODES OF INTERCONNECTED LATTICE TRUSSES 
USING CONTINUUM MODELS, PART I 


A. V. Balakrishnan 


Summary 

Continuum models are constructed for interconnected beam-like lattice trusses 
typified by the NASA-LaRC Phase Zero Evolutionary Model. For the main bus as well 
as the appendages we use equivalent one-dimensional Timoshenko beam models leaving 
only the antenna structure as lumped. The dynamic equation is cast as an abstract wave 
equation in a Hilbert space with a mass-inertia operator, a stiffness operator and a control 
operator. One novel feature is the introduction of “linkage conditions” to take care of 
interconnection of trusses. Formulas are developed for modes and mode shapes — they 
take the form of roots of determinants of matrices, albeit involving transcendental 
functions. One immediate use of the formulas involves the study of asymptotic modes. 
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1. Introduction 

This paper is Part I of a two part report and represents a continuing systematic 
attempt 11-5] to explore the use of continuum models — in contrast to the Finite Element 
Models currently universally in use — to develop feedback control laws for stability 
enhancement of structures, particularly large structures, for deployment in space. We shall 
show that for the control objective, continuum models do offer unique advantages. 

It must be admitted of course that developing continuum models for arbitrary struc- 
tures is no easy task. Attempts are beginning in this direction, nevertheless — see [6]. In 
this paper we take advantage of the special nature of current Large Space Structures — 
typified by the NASA-LaRC Evolutionary Model |9] which will be our main concern — 
which consists of interconnected orthogonal lattice trusses each with identical bays. For 
beam-like lattice trusses, an equivalent one-dimensional Timoshenko beam model has been 
developed in 171. Using this approximation, we develop an almost complete continuum 
model for the Evolutionary structure. We do this in stages, beginning only with the main 
bus as flexible and then going on to make all the appendages also flexible - except only 

for the antenna structure. 

Based on these models we proceed to develop formulas for mode frequencies and 
shapes. These are shown to be the roots of the determinant of a matrix of small 
dimension compared with mode calculations using Finite Element Models, even though 
the matrix involves transcendental functions. The formulas allow us to study asymptotic 
properties of the modes and how they evolve as we increase the number of bodies which 
are treated as flexible — as we shall see the asymptotics in fact become simpler. 

Our treatment is substantially different from extant approaches to modal analysis, 

e.g., 1 81. 

We begin in Section 2 with a brief description of the NASA-LaRC Zero Phase 
Evolutionary Model. In Section 3 we describe the one-dimensional equivalent Timoshenko 
beam model of a lattice truss, following Noor e. al. 171. In Section 4 we develop continuum 



- 2 - 


models of the Evolutionary Model in three stages: First we model only the bus as 
flexible; in the second case we model the bus as well as the laser tower as flexible; and 
finally the bus and the tower as well as the appendages are modelled as flexible with only 
the antenna as rigid. In Section 5 we develop formulas for the mode frequencies and 
shapes for all the three cases. A study of the asymptotic modes and mode shapes 
is presented in Section 6 drawing on the formulas in Section 5. The closing section, 
Section 7, contains some conclusions based on the study. 

In Part II of this paper we shall present results of numerical computations of modes 
and mode shapes based on the formulas herein; and compare them with extant calculations 
based on Finite Element Models. 
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2. The NASA-Langley “Evolutionary Model” Structure 

A schematic of the Evolutionary Model, consisting of a long truss bus and several 
appendages with varying degrees of flexibility, is shown in Figure 1. The main truss bus 
structure has 62 bays, each being a 10-inch cubical bay. The vertical appendage (Laser 
Tower) is a truss with 11 bays. There are four horizontal bay appendages each with 10 
bays (to which suspension cables are attached). There are 4 bays on the reflector tower. 
The reflector has eight 0.25-inch thick (aluminum) ribs which taper in width from 
2 inches to 1 inch over their 96-inch length. For more details see [9]. The relative 
positions of the appendges are schematized in Figure 2: s 2 , s 5 locate the horizontal 
appendages, s T denotes the tower truss; the antenna is at L; 0, s T , s 4 are co-located 



Figure 1 
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3. An Equivalent 1-D Timoshenko Beam Model of a Lattice Truss 

Here we follow Noor et al. [7] in their technique for constructing an equivalent one- 
dimensional continuum model of a Lattice Truss as an anisotropic Timoshenko beam. 

The element properties of a generic truss are shown in Figure 3. 

Let the truss axis be the x-axis and let the z-axis be the vertical, and the x-y plane be 
the horizontal, plane. Let u, v, w denote the displacement along the axes at the bay vertices. 
Let s parametrize the position along the bus axis, 0 < s < NZ, where N is the number of 
bays. Let «(s), v(.v), w(j) denote the displacement at s for the equivalent Timoshenko 
beam, and let <J>j(.y), <|> 2 (-r). be rotation angles about the x, y and z axes respectively, 
0£s< N£. Then the Timoshenko variables are related to the node displacements by: 


u(U) 

y(U) 

w(kZ) 



<M**) = ^ + v [kz, |, -|) 

- w [ike, |, =|) - W (u, =|,|) - v(*«, =|, ±1) - v(kZ, |, |) ] 

<M**) - ^ [« =|, |) + u[kt , |,|) -«(**,= |,=|) 

<d 3 («) - ^ [*«(«. f.=D + * («•!•§ -“(«•=!’ D] • 
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TABLE: Element Properties 



Battens 

Longitudinal 

Bars 

Diagonal 

Bars 

Cross Bracing 
in Battens 

Length L 

b 

l 

d 

6 

Sectional 
Area A 


A t 

A d 

^5 

Elastic 
Modulus E 

E b 

E t 

E d 


Mass Density p 

P b 

Pi 

Pd 

P5 

Element 
Mass = pAL 

m b 

m £ 

m d 

m 5 

Element 

Stiffness = EAIL 

s b 


s d 

5 8 


Figure 3 
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The anisotropic Timoshenko equations between nodes (discontinuities) are, introducing 
now the time variable t, so that 


u(s, t) 
\{s, t) 


f(s, t) 


w(j, t) 

<M*.0 

<h(S, t) 

0 


Mo 


57^ - A 2 + i4, + A 0 f(t, s) - 0 , 


a ,2 - -2 - -i ds 

where s- represent nodes, with the convention: 


s i < s < s u 1 


(3.1) 


- 0 : 

•*2 : 

5*3 ™ Sj' 

s 4 : 
s 5 : 

st> *= L \ 

A 2 = 

At = 


sensor/actuator 

appendage 

tower/sensor/actuator 

sensor/actuator 

appendage 

antenna/sensor/actuator 

C, 0 
0 C 3 


0 C 2 
-C 2 0 


M 0 


i4 0 = Diag. [ 0 , 0 , 0 , 0 , c 55 , c 44 ] 


mn 0 00 0 0 

0 m 2 2 0 0 0 0 

0 ftl33 0 0 0 


0 0 m44 0 0 

0 0 0 m 5 5 m 56 

0 0 0 /njg mgg 
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where 



c It 

C ^4 

^15 

c, = 

Cu 

C44 

c 45 


<- 15 

C45 

C55 



0 

~C\5 

^14 

c 2 - 

0 

~^45 

C 44 


0 

^55 

C 45 


c 66 

^3 6 

c 26 

c 3 = 

C36 

C 33 

c 23 


C26 

c 23 

C 2 2 


The mass coefficients m iy in the Timoshenko equation are given in terms of the bay 


parameters by: 


mn “ m 22 = m 33 


Am h + 4 m t + 4 m d + m h 

F~ 


m 4 4 « 2m$s - 2^6 


fc(8mg + 12m t + 8 m d + mg) 

6 p 2 


tm s 

mse - - 12 ^2 

The stiffness (flexibility) c iy are given by: 

4gy d li 2 

c n = 4^S t + s d + S b (l + H 2 ) 

C i 4 £21 _ llS b S d 

C 44 = — ' Cj 5 ■= - + S b (l + ) 1 2 ) 

e 3 s t 

C22 = <'33 = ^2 + 4(S d + Sg(* + ft z )) 

- *%S4 

C23 = 4(S d + S b (l + |i z )j 

&s b s d 

c 6 6 - 2 c 26 - -2c 3 6 = ^ 2 (S d + s b (t + ft 2 )) 
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where 



Evolutionary Model Parameters 

For the evolutionary model, the coefficients specialize to: 


m \\ 

= m n 

* W 33 * 1.076 

x 10 " sluglet/inch 

m 44 

* 48.31 x 10 3 sluglet-i 

inch 

m 55 

* ™66 

= 24.15 x 10 " 3 

sluglet-inch 

c \\ 

- 62.45 x 10 5 lb 


c 22 

“ C33 

= 7.06 x 10 5 lb 


C44 

= 353. 

14 x 10 5 lb-inch 2 


c 55 

“ c 66 

- 1540.46 x 10 s 

lb-inch 2 . 
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4. Continuum Models of the Evolutionary Structure 

We develop now (flexible) continuum models of the evolutionary structures at levels 

of increasing complexity: 

i) Bus only as flexible 

ii) Bus and tower as flexible 

iii) All (bus, tower and appendages) as flexible, only reflector lumped. 

In all cases we shall obtain the generic model dynamics as an abstract wave equation 

in a Hilbert space: 

Mx(t) + Ax(r) + Bu(t) = 0 

where 

*(•) e Hilbert Space H. 

M is the mass-inertia operator: At is a self-adjoint and positive definite 
linear bounded operator on X onto K with bounded inverse 

A is the stiffness operator: closed-linear operator with domain dense in X: 
self-adjoint and nonnegative definite with compact resolvent 

B is the control operator: B maps finite-dimensional Euclidean space into H 

u( ) denotes the control (input) . 

See 11] for the first development of such a model. Among the advantages of this generic 
formulation is the close similarity of FEM and truncated modal models - excepting only 
for dimension not necessarily finite! We begin with the first case: 

Case 1: Bus Only as Flexible 

In this model the tower, the appendages and the reflector are modelled as offset 
lumped masses, as are the controllers, and the bus represented by the equivalent 1-D 
anisotropic Timoshenko model. Let denote the location of the lumped masses. It is 
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convcnient at this point to invoke the abstract or function space representation as in 
[1]. Our function space denoted H is taken as: 

X - Z^[0, L} 6 x /? 6x6 

with elements denoted x : 



«(*) 

v(j) 

w(s) 

f(s) « , 0 < s < L 

<M*) 

<M*) 

<M*) 


/( 0 ) 
f(s 2 ) 

/( J r) 

b - 

/fo) 

m 

and the norm in X is given by: 

L 

IMI 2 - J ll/OOII 2 ds + ||6|p . 

0 

We now define the operator A : 

A f - S 

b c 

where 


(4.1) 


(4.2) 


g(s) - -A 2 f"(s) + A\ f\s) + i4 0 /(s). 


^ s k S U j » / ** . 
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-L,/(0) - A 2 f\ 0) 

A 2 {f\s 2 -) ~ /'(*a+)) 

Mf( s r) - /'(■*,+)) 


Z.i/(Z.) + A 2 f{L) 


( 4 . 3 ) 


where 




0 — C 2 
0 0 


The domain of A consists of functions which are continuous and piecewise smooth: in 
fact are in 7t\s„ s 1+1 ), i - 1. •••, 5, and the first derivative is possibly discontinuous 
at s = $;. A is then a closed linear operator with domain dense in X and is self-adjoint 
and nonnegative definite. Moreover for x in D(A) 


[Ax, x] 



0 


L 

r 

* 

«'(*) 

| 

c, 

v'(s) - <h(s) 

0 


w\s) + <t> 2 ( 5 ) 


ds 


u\s) 

v'(j) - <t>3 (s) 

W'(S) + (jljfa) 


■ 

L 

* 



<t >i(s) 


ds + 


c 3 

<H(s) 

i 

<t>2 (*) 


C 

) 

„ 

4 > 3 ( s ) 


<J>3(5) 

- 


ds 


( 4 . 4 ) 


( 4 . 5 ) 


where 

0 - 1 C 2 

0 0 

and the potential energy of the beam 

[Ax , x\_ 

2 

It is of course assumed that C, and C 3 are positive definite and nonsingular, and hence 
A is nonnegative definite. 


C^ 0 
0 C 3 
0 0 
-c\ 0 
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We have thus obtained our “stiffness” operator. Next we need to define the mass/ 
moment operator M. 

f Mof 

M 

b M b b 

We proceed now to define M b . M b is “diagonal”: 

M„b i 

M b b - (4.7) 

6 

where 

bi 

b = 6, g R 6 

b 6 

and M b . are nonsingular, symmetric and positive definite. 

Finally we define the control operator B . Figure 4 is a schematic of the Evolutionary 
structure showing the disposition of the force actuators and the corresponding axes along 
which they act. There are 8 actuators. Hence let U denote the 8x1 column vector: 



Then 


where 


«i 
« 2 

«8 

0 

Bu = x ; x = 

B V U 

b i 
b 2 

B U U = = b 

b 6 
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Figure 4 

Schematic of Evolutionary Structure 
Showing Disposition of Actuators/Sensors 
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«3 
"4 
0 

62-O - 

R t x 

b$ = 0 “ 

*L 

where, with r denoting the position vector, r(s) denoting position vector along bus axis: 

Ro = ^controller) - r(0) 

R t = r(controller on tower) - r(s T ) 

R, = ^controller at 1 = j 4 ) - rfo) 

R l - r(controller at s = L) - r(L) . 

We note that the numerical values are: 

R 0 - 0 

0 

R t = 0 

100 

/?, - 0 

0 

r l = 0 . 

40 
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Hcnce 


B*x * Bjj b 


v(0) 

M<0) 

u(s T ) + 100<|> 2 (5 T ) 
v(s T ) - 100<t>,(s r ) 
v(j 4 ) 

M<5 4 ) 

u(L) + 404> 2 (L) 
v(L) - 40<J»j (L) 


Case 2: Bus and Tower as Flexible 

For this case and the next it is convenient to change notation slightly. We use 
f(x, y, z) in place of f(s), so that 

As) - As, 0, 0 ) 

denotes the displacement vector along the axis of the bus and 

f T (s) = /(J T , 0, s) , 0 <s<L r 

will denote the displacement vector along the axis of the tower truss in the equivalent 
1-D Timoshenko model, with L r denoting the length of the tower. Since the tower truss 
axis is now the z-axis, we redefine the tower truss coefficient matrices using a subscript: 


’2.7 


'1,7 

0 


0 


'3.7 


■"2,7 


M ,7 


-C 


2,7 0 


'0,7 


= Diag. [ 0 , 0 , 0 , c 44 , c 55 , 0 1 
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ftl 33 0 

0 

0 

0 

0 


0 

m 22 

0 

0 

0 

0 

Mq T “ 

0 

0 

m n 0 

0 

0 


0 

0 

0 

"»66 

ms6 

0 


0 

0 

0 

™ 56 

"»35 

0 


0 

0 

0 

0 

0 

"144 

L, T . 

0 

0 

-C 2 ,r 

0 






c 55 

C45 

c 15 




Jf} 

H 

II 

c 45 

C44 

c 14 




i 

c 15 

c 14 

Cll 




i 

C45 

-<55 

0 




N 

fc. 

C44 

"^45 

0 




i 

c 14 

“^15 

0 





C22 

^23 

^26 




Q.r “ 

c 23 

c 33 

c 36 

• 




C26 

c 36 

c 66 





In the abstract version, the Hilbert Space X now is given by 

X - LjIO.L ] 6 x LjfO ,Lj.] 6 x R 6x6 


x 


f(s, 0,0), 0<s<L 
f(s T ,0,s), 0 <s<Lp 

b 


where 

/( o, 0, 0) 

f(S2, 0 , 0 ) 

f( s T’ ^ t ) 
/(*4, 0, 0) 
/(*. o, 0) 

f(L, 0 , 0 ) 
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IWI 2 - J i/(*. °. °)i 2 * 
0 


+ J T \ns T ,o,s)\ 2 ds + in 2 . 
0 


The domain of A consists of functions 


/(*,(), 0), 0 <s<L 
f(sj ,0,s), 0<s<Lf 


where 

f(s, 0, 0) and f'(s, 0, 0) 

are absolutely continuous and f\s, 0, 0) has an ^-derivative in the sub-intervals 
0 < * < s 2 , s 2 <s <s A , s 4 <s<s 5 , s 5 <s<L ; 


and 


f(s T , 0, s) and fX s T> 0* J ) 

are absolutely continuous with f\s T , 0, s) having an Lj-derivative in 0 < s < Lj 
Moreover the following “linkage conditions” are satisfied: 


i) f(s, 0, 0)|,. Xj- •= /(*r» 0, j)|,.o 

ii) L t T f(s T ,0,0) - A 2T f z (s T ,0,0) + A 2 (f x (s T -,' 0,0) - f x (s T +,0,0) « 0. 


The stiffness operator A is now defined by 

/( •, 0 , 0 ) 

x = I f{s T , 0, ) 

b 


Ax 


g('.0, 0) 

g( s T’ ’) 

V 


where 

0, 0) - -A 2 f IX (s, 0, 0) + A x f x (s,0,0) + A 0 f(s, 0,0) 

0 < s < s 2 , s 2 < s <s 4 , s 4 < s < s s , s 5 < s < s 6 , 


g(s r , 0, s ) 


~A 2 f f zt (s T , 0, s) 


+ Aj T f t (s T , 0, s) + A 0 j f(s T , 0, s) 


0 < s < Lp , 
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-L,/(0, 0, 0) - A 2 f x { 0, 0, 0) 

MfM- 0 , 0 ) - f x (s 2 +, o , 0 )) 

A f j ,rf( s T> + ^2jfi^ s T' ^t) 

* 2 (/,(*-, 0, 0) - /,(s 4 +, 0, 0)) 
o, 0) - f z (s 5+ , 0 ,0)) 

Lif(L, 0,0) + A 2 /,(L,0,0) 

Thus defined, it is easy to verify that A is closed, self-adjoint and nonnegative definite 
and that 

jcl 

J — 2 — L “ Elastic Energy of Bus + Elastic Energy of Tower Truss . 

Finally the mass/moment operator M is defined by 


Mx 


Afo /(•,(>, 0) 

Mq j f(s f, 0 , •) 


M b b 


where 


M b b 


M bi b i 


and M b are positive-definite and nonsingular. 

Finally we define B the control operator. First we define 

0 


BU 


x = 


B„V 


B U U 


col.[0, u u u 2 , 0, 0, 0] 
0 

col.[u 3 , u 4 , 0, 0, 0, 0] 
col.[0, u 5 , u 6 , 0, 0, 0] 
0 



• 

u 7 

- 

col. 

u 7 , Ug, 0, R l x 

«g 



■ 

0 

« 
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B*U - B*b 


v(0, 0, 0) 
w(0, 0, 0) 
u{sj, 0 , Lij-) 

v(s r , 0, Lj) 

\{s 4 , 0, 0) 

wfo. 0, 0) 

u(L,0,0) + 404> 2 (L,0,0) 
v(L,0,0)- 404*! (L, 0,0) 


Case 3: Bus, Tower and Appendages Flexible 

We now generalize to the case where the main bus, the laser tower and the hori- 
zontal appendages are modelled as flexible lattice trusses — or more precisely, their 1-D 
Timoshenko beam equivalents. We shall be briefer in our descriptions since we will 
follow the pattern already set in Case 2. 

Thus let the subscript s 2 denote the coefficient matrices for the appendages at s = $2 
and similarly the subscript *5 for the appendage at s = s 5 . Then 

Cl* 0 

' 0 c„, 

I 0 C «. I 


A o s - Diag. [ 0 , 0 , 0 , c 55 , 0 , c 44 ] 


L 


l.s, 


0 ~^2,5 1 

0 0 


m o,s 2 


m 22 0 0 0 0 0 

0 m\\ 0 0 0 0 

0 0 m33 0 0 0 

0 0 0 m 5 5 0 m 56 

0 0 0 0 m 44 0 

0 0 0 ms6 0 "*66 
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c 44 Cj 4 C 4 j 

^\,Sj “ c 14 Cll C 15 

c 45 c 15 C 5S 

-C45 0 C44 

"<15 0 C14 

~<55 0 C45 

c 33 ^36 C 23 

Q,S 2 “ c 36 ^66 ^26 

c 23 c 26 C22 

A 2,S , “ ^2,5j 

•^1.5, “ ^l.Sj 

^O.S, “ ^O.Sj 

^0 .S, “ ^O.Sj 

^ 1 , 5 , “ ,S 2 • 

The appendage displacement vectors are then 

f(s 2 , s, 0) , -ti < s < ti 

f(s 5 , s, 0) , -t 2 < s < t 2 

for the evolutionary truss £1 = £ 2 • Thus let 

f(s2,s,0), -£i<s<£, 
f(s T ,s, 0), 0<s<L r 
f(s s ,s,0), -£ 2 <j <£ 2 
f(s, 0,0), 0 <s<L 
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/(O, 0, 0) 
/(^2,-ti.O) 
f(s 2 , +ti, 0) 
f(Sjy 0, Lff) 

f(s A , 0, 0) 
Ms. -l 2 . 0) 
Ms, +t 2 , 0) 
f(L, 0, 0) 


-L,/( 0, 0, 0) - A 2 f x (0, 0, 0) 

-L x Si f(s 2 , -t\, 0) - A 2 Si f y (s 2 ,-t\,0) 
L| ,s,/( 5 2» +^1» 0) + ^2 ,S 2 fy( s 2' 

Lj T f(,s T , 0, Lj) + A 2 jf z (s T , 0, Lj-) 
A 2 (/ x (5 4 - 0, 0) - f x (s 4 +, 0, 0)) 

-L x Si f(Si, -Z 2 , 0) - A 2 Si f y {ss,-t 2 ,0) 
L x 'S t f( s S> +^2 » 0) + A 2 s Jy(S5, +Z 2 , 0) 
L,/(L, 0, 0) + A 2 f x (L, 0, 0) 


Plus Linkage Conditions: 

( 1 ) f(s, 0, 0)|,., r - f(s T , 0,s)|,.o 

-L, r f(s T , 0, 0) - A 2 T f z (s T , 0,0) + A 2 {f x {s T - 0,0) - f x (s T +,0,0)) = 0 

(2) f(s, 0,0)|,. fl = /(S 2 .^ 0 )|,. 0 

^2(/,(^2-,0,0) -/,te+.0.0)) + *2.s a (/yta.<K 0) - fy{s 2 ,0+,0)) = 0 

(3) MO, 0)1,.,, = Ms.*, 0)|,. 0 

A 2 {f x (s 5 -, 0,0) - f x (s s +,0,0)) + ^ 2> s 3 (/y( J 5,0-,0) -/ y (^,0+,0)) = 0. 
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Remark 

If suspension “ends” are treated as “free-free,” then remove 

/C*2, -*i, 0) 
f(s 2, +t|, 0) 

from 6 and instead take: 


^l,S,f( s 2< ~f-l> 0) + ^2,S t fy( s 2t -ti, 0) - 0 
+^1. 0) + ^2,S t fy( s 2> + ^1» 0) “ 0 
and similarly for the other suspension beam. 

Finally: 


x - 


/ 

b 


X - m-U, t, ] 6 x mO.Lr ) 6 x LjMj, t 2 ] 6 x LjIO.L] 6 x tf 6x8 


where 


Ax * y 



y « 


s 

c 


S(*2. *, 0) 

Sfar* °> ) 

g(S5, ,0) 
g(‘, 0,0) 


£(- v 2> s > 0) “ A2 s i fyy(S2,S,0) + i4j s > 0) + A) ,S t f(. s 2r S, 0) 

-£i < s < tj , 

g(Sj-,S,0) — ~^2 ,T fn^T ' S ' 0) + "^1 j fi( s T> S ' + ^0,Tf(- s T> 0) 

0 < S < L , 

g(s 5 ,s,0) = -4 2 5j / yy (j 5 , j, 0) + A l Si fy(s$, s, 0) + i4 0 Sj /(j 5 , 5, 0) 


-t 2 < J < t 2 , 
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g{s, 0, 0) - -A 2 f xx (s, 0, 0) + AJ x {s, 0,0) + A 0 /(i,0,0) 

c - A b f ; / subject to linkage conditions . 


Then A is self-adjoint and nonnegative definite and 

id|i *1 „ [Sum of Elastic Energy of Tower, Suspensions and Main Beam] 

Next, the mass/inertia operator M is defined by 


Mo'SjfiSi* '* 0 ) 

Mx ■* I ■, 0) 

M 0 /(-,0, 0) 

M b b 


where again 


M„B 


M b b , 


where are mass/inertia matrices. 

Finally we define the control operator B. 

I / 

x - 

I b 


BU 


0 

Byf/ 


By U = col. I B, , b 2 , B 3 , B 4 , B 5 , b 6 , bn , Bg ] 
Bj = col. 10, «i , «2 » 0 , 0 , 0 ] 


B 2 = 0 - B3 

B 4 - col. [ «3 , u 4 , 0 , 0,0,0] 
b$ = col. [ 0 , «5 , «6 . 0 , 0 , 0 ] 


B 6 = 0 - bn 
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b* 


«7 

"8 

0 


R, x 


« 7 
w 8 
0 


Hence: 


B*x - BJ6 - 


v(0, 0, 0) 

MO, 0, 0) 
u(s T , 0, Lj) 
v(sj, 0, Lj-) 
v(j 4 , 0, 0) 

H*4. 0, 0) 

u(L, 0,0) + 40<t> 2 (L,0,0) 
v(L,0,0) - 404>,(L,0,0) 
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file: june91y 


5. Mode Formulas 

In this section we develop formulas for modes and mode shapes to find the modes 
we need to solve the eigenvalue problem: 


Ax - (o 2 Mx . 


(5.1) 


Letting 



we begin with Case 1. 


Case 1: Bus Only Flexible 

In this case (5.1) translates into: 


x 


f 

b 


-A 2 f"(s) + A 1 /'(s)+ A 0 f(s) - CD 2 M 0 f(s) , 

A b f - a ?M b b 


*, <•*< s u l * 


where the second equation can be expanded as: 

-L,/(0) - A 2 f\ 0) - M m /(0) 

- M bti m, i « 2 5 • 

f - l /(^) + ^ 2/(^0 “ ^b , L m j 


For combining (5.2) and (5.3a) let 


j4(g>) - 


0 

A 2 (Ao-(& 2 Mo) 


I 

A 2 Ai 


(5.2) 

(5.3) 


(5.3a) 


(5.4) 


(a 12x12 matrix). Let 
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where, of course, 


Let 


Then (4.9a) yields: 




PM P n {s) 

Pn(s) P 22 (s) 


Pii(s) - Phis) 
Puis) - Puis) . 


Zi - A 2 ((a 2 M b Q - L\) 

- <&AiM b i , i = 2, .... 5 

^6 “ A 2 {(0 2 M b l - L\) . 


(5.5) 


M+i ) 




/(*,) 

fiSi-)-*ifiSi) 


with the convention that 


/'(*!-) - /'(0-) - 0 


and condition (4.9a) requires that 


/'(L) - ^/(L). 


But we can write 


where 


f\L) - Ztf(L) = D(to)/(0) 


D(to) « ^(L, -wX.l) 


~<o i A 2 M b 'S 


0 

/6 


a(co)(i,-* 4 ) 


f 6 0 

-W 2 A 2 M b 4 I b 




~® 2 AiM bl 


d(ia)(s- i -s 1 ) 


/ 6 

-01 1 A 2 Mi 


b.2 


0 
f 6 


a(m)»2 


(5.6) 


0 

/6 


^2(^1 — v?M b ' 0) 


(5.7) 
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where / 6 is the 6x6 Identity matrix. Thus the mode frequencies are determined from 


P(G»bet = 0 

and the mode shapes from the corresponding eigenvector /( 0): 


O(co)/(0) = 0 

the corresponding /(*,■) being determined from (5.6). Or, more explicitly 
f(s) = |f 6 0 | e AimUs ~ s i ) T { • r,_. 


(5.9) 




/( 0 ) 

Ali-Li -to 2 M M )/( 0) 


s < 5, 


(5.6a) 


where 


U 

| -0) 


0 

* 

/6 


i - 2,3,4,5, A;. - s, + i - ■r, • 


We have thus “reduced” a mode determination problem to finding the zeros of a transcen 


dental function 


P«*»l D e. “ °‘ 


The crucial calculation is that of the matrix exponential We note that we 

can “expand” D(co) as: 

6 2 * 

D(g>) = I to D k 


since 


D(0) = 0 . 



-30- 


Pure Modes 

The evolutionary model trusses are actually isometric: 

c ij - 0. i*j 

m ij - 0 , i * j 

M bo - M bt - Diag. ( .05 , .05 , .05 , 1.9 , 0.95 , 0.95 ) 

0.28 0 0 0 0.71 0 

0 0.28 0 -0.71 0 0 

Af tj - - 0 0 0 0 0 0 

0 0.71 0 1538 0 0 

-0.71 0 0 0 53.9 0 

0 0 0 0 0 1494 

0.18 0 0 0 13.23 0 

1 0 0.18 0 -13.23 0 0 

M, t - 0 0 0.18 0 0 0 

0 13.23 0 1132 0 0 

-13.23 0 0 0 1132 0 

0 0 0 0 0 7.3 

0.38 0 0 0 22 0 

0 0.38 0 -22 0 -0.91 

M bf “ m l - 0 0 0.38 0 0.91 0 

0 22 0 1511 0 120 

-22 0 -0.91 0 1459 0 

0 0.91 0 120 0 229 

Mo - Diag. ( 1.08 xl0'\ 1.08xl0 -3 , 1.08xl0" 3 , 48.3xl0" 3 , 24.15xl0 -3 , 24.t5xl0 -3 ) 

A 2 - Diag. (62.45x10 s , 7.06x10 s , 7.06x10 s , 353.1x10 s , 1540x10 s , 1540x10 s ). 

The inertia matrices M b . are nearly diagonal. If we retain only the diagonal terms, 
we can easily see that there are “pure” modes: a pure “axial” mode in which 
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f(s) - a(i) 


1 

0 

0 

0 

0 

0 


and pure “torsion” mode: 


f(s) - a(s) 


and we can calculate the corresponding mode frequencies (and shapes). Thus let 


Then we have 


D(co)c i = d( o))ej 


where 


d(v>) 


9 


Pn(As) Pn(*s) 


-0? m b,L . 

. 


• 

c 11 


P 2 i(A s ) /Ws) 



c i j 


/ > n(A i ) 

/ > 12(A | ) 


1 

0 

P2t(A ; ) 

/ > 22(A 1 ) 


C 1 1 

l 

f’nCAi) 

P>2(Ai) 


1 



P 2 2(Al) 


-<0 2 m b , 
c 1 1 



0 

1 



where 


where 


and 


and the mode shape 


a(s) - 1 0 | • 


- s ,+ 1 _ J i . i m 1, 

P 1 1 CO — cos Aj 
Pn(s) - 

/*2i(0 - -A, sin Xs 
P 22(0 - cos A s 

X - (oVmi |/c] i 


- 

the 1-1 entry in 

M b ,i 

P\ l(S-Si) 

Pn(s-Si) 


1 

Pl\(S-Si) 

Pn(s-Si) 

• 



ct i 


Pn( A,) 

P 12( A l) 


1 

^i(Ai) 



-a?m bl 

Cl 1 


We list below the first few modes corresponding to 



-33- 


Pure Axial Modes (Hz) 
29.1 
82.45 
116.35 
187.4 
218.7 

281.9 - 
370.3 

610.9 


Pure Torsion Modes 
Here 


where 


f(s) - a(s)e 4 


e 4 


0 

0 

0 

1 

0 

0 


D(( 0)e 4 - ds(<o)e 4 


d(a) - 




Pn(A 5 ) P 12 (^ 5 ) 


«6.4.t j 

. 


• 

C 4 4 


P 2 i(A 5 ) P 22 (A 5 ) 



1 

4,5 

r 4 4 


0 

1 


Pi. (A,) 

P.2(A ( ) 


1 

0 

p 2 .(a,) 

P 2 2(Ai) 


-® 2m fc.4,i 
C 4 4 

1 

P..(A,) 

Pu(Ai) 


1 


P 2 i(Ai) 

P 2 2(Al) 


-® 2 m fc( 4.1 



C 44 
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whcre 


m b, 4 ,i m 4x4 terms in M b i 
X - coV^7 
P\\ (j) - cos Xr 
fn(s) - ^ 


/ > 2 i(r) = -X sin Xs 


P 2l( s ) m COS Xs 


The first few modes are: 


Pure Torsion Modes (Hz) 
1.2 
4.29 
6.79 

30.6 

41.6 
67.46 
94.53 

127.4 

163 

186.9 

208 

232 

249 

292 

317 


Case 2: Bus and Tower Flexible 
For this case (5.1) yields: 


D(o>)/ - 0 
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where 


/ “ 


D{ ©) 


f z ( s T * S ) I j-0 

/( 0 , 0 , 0 ) 


Du D 12 
D 2 1 D 22 


12x1 


O ,, ■= | a 2,j ( l i,t ~ 1 e T / 


0)2 = ,T ~ I e T 0 


I 0 I ji(tD)(Sj.-Sj)- 


i«o)Si 


^(-Lz -W 2 Wfc f i) 


Djl , | ^(L, - / | ■ (e t ' m) ' L - S ' ) T s ) ■ 


i ( to )( S 4 - 5 T ) 


0 

2^2.T 


A^A 


D2! . | fa 'I ■ ^ wtL - s ' ) T s ) ■ (« J, “ ),S ‘- S<, n) 


J(oi)(S 4 -S r ) 


/ 0 
^2 0) ,r ^ 


(c J '( a>KS r _S2> T , 2)e' <<a>>52 


A?.(L, 


where 


lO((0)| 


^(CO) 


T: 


| 0 |l 022 ~ O11O21O11D22 


0 / 

I A2 (A 0 - w 2 M 0 ) A2 Ai 

/ 0 
-V?AiM b<i I 


- 2, 3, 4, 5 


I 

-a >X. 
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J T ((0) 


0 


/ 


A 2,t( A 0,T~ rfMoj) A 2,T A \,T 


V - 


/, 

fo 


D((o)f - 0 


Mode Shapes: Tower 

f(s T , 0 , z) - |/ 0 |e' r(< “ )r 


X 


I 


o 

f t (s r , 0,0) 


/ 0 
0 0 


e J(a»A Ir ^ ^(a»A, 
0S2<Ly. . 


/(0, 0, 0) 

^(-L,-(oX,i)/(0,0,0) 


Case 3: Bus, Tower and Appendages Flexible 
For this case (5.1) reduces to: 


£>(«>)/ = 0 


where 


/ * 


fyMi, o, 0) 
/ y+ te. 0, 0) 

f z ( s T> 0 * 0 ) 


fy-(S 5 , 0 , 0 ) 

/y + C*5. 0, 0) 
/(0, 0, 0) 


36x1 


D(g>) - {£>■;} 1,7-1 6 
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£>n - | 

-L 1>Si -(0 2 M Sj 


£>22 “ 



£>33 “ 

L\ j-G> 2 Ms 2t 


£>44 = 

-L, Sj -G) 2 A/ s 


£>55 = 

| £l .5,-^5, 

£>16 = 

| -L, Sj -(0 2 / 

-*2.5, 

£>26 = 

|l, iSi -(0 2/ 

*2.5, 


0 

/ 


0 

/ 


/ 0 
0 0 


)A, 


—Ali (L\ + (S?M j, o ) 


/ 

0 1 i(w)A, 

0 

0 1 


-A2 (L\ + (0 2 Af b o ) 


£> 36 


Lj j - <o 2 / 4 2 ,r « J( “ )rLr ‘ 


/ 0 
0 0 


i(«)(A,+ A,) 


-v 4 2 (£i+g> 2 M m ) 


^>46= |(-£l.S } -0>X,.-« 1 ) -^2.S s | eJ,<t0)i5( ‘ tj) • 


£ 0 J(»)A 4 ~ i(<a)(Aj + Aj + A,) 

£ i 4 e 

o o 


-Aj (L\ + w 2 M fc o) 


£>56 = |(£M 3 - ft)2W S s .* I ) A 2,s 5 |^ ( “ )sjtj ' 


/ 0 
0 0 


^l(c#)A 4 j, ^«0)(A, + A, + A,) 


-A 2 (£i +o» 2 Mfc,o) 
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066 


U, -v?M b , A |^(«)< a » 4 a 4) t ^(«)A, . 

/ 0 

1 * | 

~ A 2 j / 


(w)(A 2 + A|) 


I 

-Ai(]Lx+ii?M bA ) 


065 


Li-afMhL 


A 2 e J( « )A ’ 


0 

~ A 2 A 2 ,S s 


06 A - j Li-<o 2 M b L A 2 | e i(c#)A * 


0 

A 2 A 2,S s 


063 - L]-(ii 7 M btL A 2 ^ ( " )<a * + A 0 r 4 ^ ( “)a» . 


~ A 2 A 2 ,T 


0, 


62 


L,-o) 2 M fc L /l 2 e' <( “ )(A>+A 4 ) r 4 *' (a,)Aj . 


-v4 2 Lj r 


0 

/ 


>(•)*» 


/ 

"^2^2.5, 


'61 


L,-(oX t M 2 e '<“»< A ’* A «> r 4 ^ (<0)A ‘ • 


A (<o)A 2 


I -i4 2 L, 7 


/l 2 /l 


2.5 a 


|D(tO)| - 0„ • 0 22 ' • • 055 • (066 ~ 0610T|016 - - - D 65 Ds\d 56 ) \ 
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6. Asymptotic Modes 

In this section we use our mode formulas of Section 5 to study the asymptotic 
behavior of modes for all the three cases. 


Case 1: Bus Only Flexible 

For the asymptotic study we use the expansion 

6 

£>(«) - I ai Lk D k . 

\ 


For large to therefore the roots of 

|D(to)| - 0 


are those of 


Pel - 0 


with increasing accuracy. Now 

£> 6 = A~2 X M b L Pi 2 (A 5 ) Al'M b ' 5 / > n(A 4 ) *" P\i(A\) A 2 X M bX 


where 


Pn(s) -| / 0 | e 


A(t»)s 


I 

0 


Hence 


m - o 


if and only if 


|/>, 2 (A-)| ■= 0 , for some i - 1, 2, 3, 4, 5 . 


But these roots are recognized as the modes of a “clamped-clamped beam. Thus 
asymptotically all the “clamped-clamped” modes correspond to every beam segment between 


nodes. In fact for large <0, 

Pi 2(A) 


6 sin X k (to)A 


where e * are the unit vectors, k — 1, ...» 6, 
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X A (co) - toVy^ 


A 2 M 0 e k - y*e t 


The eigenvalues of / > i 2 (A i ) are given by 


sin Vft wA ( 


A - 1, 6 


or, the modes are given by 


(OA, - nK 


® " a i — t r * 1 , . . . , S , ^ 

A| Vy* 


For the evolutionary model 


w “ 2T^ c kk /m kk 


V “ 2A ~^k' m kk Hz 


For A 4 the largest segment, this yields for the axial mode: 


and for the torsion mode 


v - (165.4)n Hz 


v - (58.7)n (A 4 ) 


(77 . 1 )n (A 3 ) . 


It is difficult to recognize these in the few modes we have calculated. 


Case 2: Bus and Tower Flexible 


Here |D(( 0 )| - 0 for large co yields 


|D,,«o)||D 22 ((o)| - 0 


|D,, (0))| -= 0 
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means the roots of 

|/ 0 \e A T ( * )L T |'| - 0 


which are recognized as the “clamped-clamped” modes of the tower truss — as we 
should expect. 


|D 22 ((0)| - 0 


yields 

j A 2 * M b i / > i 2 (A4) i4 2 M b 4 | = 0 

= 0 


The first relation is equivalent to: 


|/» i 2 (A 5 )| = 0 ~ clamped-clamped modes of segment A 5 

| J P l2 (A 4 )| = 0 ~ clamped-clamped modes of segment A4 

The second relation yields 

|/» 12 (A,)| - 0 ~ clamped-clamped modes of segment Aj 


and 

|(f > n(A 3 )-/ , i2(A3)Ai 1 A i r )Pi 2 (A 2 ) + Pi 2 (A 3 )P 22 (A 2 ) | - 0 


Since asymptotically 


6 

P\ i(A) - Z cos X k (( 0 )Ae k 
1 

6 

P 22 (A) * Z COS X k ((0)Ae k . 
1 


Now 

A 2 j€ k = 0 t k * 1, 2, 3, 6 

(corresponding to “displacement” modes about the tower axis) and the torsion mode and 
hence for k *= 1 , 2, 3, 6 
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(^ 11 (^ 3 ) - Pnib-fiAi' j)P\i(bi)e k + Pn(&i)Pn{hi)e k 


ff , sinX t «o)A 2 ^ sin X k «o)A 3 „ , w 1 

- I [«» **(»)* — x^j-J + — x^T' cosX * ( 0))A2 J <? * 


sin X t ((D)(A 3 +A 2 ) 

X^j e *- 

Hence we see asymptotically the clamped -clamped displacement modes of the segment 
(A 3 + A 2 ). Hence we have the clamped-clamped displacement modes of segments: 

As 

A 4 

A, 

A 2 + A 3 

but these are now recognized as the segments between lumped masses. We note that for 
the evolutionary truss 

A 2 + A 3 *= 205 < A 4 . 

Hence these modes are still too high. 


Case 3: Bus, Tower and Appendges Flexible 
Here 

|D(( 0 )| = 0 

asymptotically 


P,, ((0)1 |D 22 (g»| |D 3 3(co)| |D 44 ((0)| |D 55 ((o)| |D 66 «d)| = 0 


H>u« o)| - 0 
P22«o)l - 0 

|D 33 (( 0 )| = 0 

|D 44 ((0)| - 0 
Pss((0)l = 0 


= clamped-clamped modes of each appendage (length £,) at s 2 
clamped-clamped modes of tower truss 
= clamped-clamped modes of each appendage at s$ 


|D 66 (<o)| = 0 => bus modes 
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- (M bfL Pn{*4 + *s)A~i' M b<A ) 

• ((Pn(A 3 )-/ > 1 2(A 3 )i4i , L 1 j-J/’iztAz+Ai) + Pu(^t)Pu(^i+^\) “ 0 
■ |Pj 2 (A 4 + A 5 )| - 0 clamped-clamped modes of segment (A< + A 5 ) . 

As in Case 2, for e k = 1, 2, 3, the displacement modes, we have 

Aj T^k “ ® 

and hence we obtain 

sin X, a (CD)(Ai + A 2 + A 3 ) - 0 . 

Or, we have the clamped-clamped modes of the segment (Ai + A 2 + A 3 ). But 

A 4 + a 5 

Aj + A 2 + A 3 

are now the segments between lumped masses. Moreover for the evolutionary truss: 

A 4 + A s = 295 
A] + A 2 + A 3 « 330 . 

The clamped-clamped mode frequencies corresponding to the segment Ai + A 2 + A 3 are 
given by 

v = -ttk Vcn/tfiii Hz ~ (Axial) - (20.5)n Hz 
oou 

and corresponding to the segment A 4 + A 5 : 

V - Vc n /m u « (22.9)« Hz 

which are now low enough to be found in the range of modes of practical interest! 
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7. Conclusions 

It is feasible to construct continuum models of flexible multibodies if they take the 
form of large interconnected trusses with many bays where advantage can be taken of 
1-D equivalent Timoshenko beam models. Using these models it is possible to construct 
formulas for modes where the matrix size is insignificant compared to the Finite Element 
version. However transcendental functions are involved. It is possible to make explicit 
use of the mode formulas to estimate asymptotic modes. Asymptotic modes would appear 
to be more realistic as the number of flexible parts which are modelled as continua 
increases. The asymptotic modes then are recognized as the clamped-clamped modes of 
beam segments between lumped masses. 
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